Three-body problem in a two-dimensional Fermi gas 
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We investigate the three-body properties of two identical f fermions and one distinguishable 
\. atom interacting in a strongly confined two-dimensional geometry. We compute exactly the 
atom-dimer scattering properties and the three-body recombination rate as a function of collision 
energy and mass ratio m-f/mj.. We find that the recombination rate for fermions is strongly energy 
dependent, similarly to what was found for two-dimensional identical bosons. For m-f < mj., the 
s-wave atom-dimer scattering below threshold is completely described by the scattering length. 
Furthermore, we examine the tti bound states (trimers) appearing at large m^/m^ and find that 
the energy spectrum for the deepest bound trimers resembles that of a hydrogen atom confined to 
two dimensions. 



In recent years, ultracold atomic gases have proven to 
be an extremely versatile system, providing elegant re- 
alizations of a range of quantum many-body phenomena 
such as the BCS-BEC crossover and the Mott transi- 
tion [1, 2]. In particular, the cold-atom system involves 
short-range interactions that can be tuned to have large 
scattering lengths, thus rendering the low-energy physics 
insensitive to the details of the interaction potentials and 
essentially "universal". This has enabled the study of 
universal few-body physics, which in turn has had major 
consequences for the many-body system. For instance, 
the scattering length of diatomic molecules (dimers) was 
necessary for a complete description of the BCS-BEC 
crossover [3-5]. 

Few-body inelastic processes in the cold-atom system 
are also important since they limit the lifetime of the gas 
and constrain the densities that can be achieved in ex- 
periment. Furthermore, they can act as an indirect probe 
of the quantum system, e.g., the first experimental evi- 
dence for the Efimov effect [6] was deduced from three- 
body losses [7]. There is even the prospect of generating 
strongly correlated phases using few-body loss processes: 
dimer-dimer dissipation has already been shown experi- 
mentally to induce correlations [8], while it has recently 
been proposed that three-body dissipation can be used 
to engineer the Pfaffian state in two dimensions (2D) [9]. 

Recently, fermionic atoms were successfully confined 
to a quasi-2D geometry, allowing the experimental inves- 
tigation of 2D Fermi gases [10-16]. However, despite the 
importance of universal few-body physics for correlated 
cold-atom systems, a complete description of the three- 
body properties of a 2D Fermi gas is still lacking. In par- 
ticular, for two identical j" fermions and one \. particle, 
there is currently no calculation for the rate of recombi- 
nation of three atoms into an atom and a dimer, while 
the existence of three-body bound states (non-Efimov 
trimers) for mass ratios m-|-/mj. > 3.34 has only very 
recently been elucidated [17]. Moreover, one cannot triv- 
ially extrapolate from the scattering problem in 3D since 
the few-body scattering properties will be strongly mod- 



ified in 2D: for identical bosons, the three-body recom- 
bination rate vanishes at low energies [18], in contrast to 
the 3D case. 

In this paper, we address this gap in the literature 
and investigate the three-body problem in a 2D Fermi 
gas. Here we compute exactly the atom-dimer scattering 
properties and the three-body recombination rate as a 
function of energy and mass ratio m-f/m^. The recom- 
bination rate has a similar dependence on energy to the 
case of identical bosons [18] , but surprisingly we find that 
higher partial wave contributions can be significant at low 
energies for both bosons and fermions. This result could 
potentially be important for any 2D experiment seeking 
to realize itinerant ferromagnetism [19]. Furthermore, 
bound states can impact the scattering properties, and in 
particular the presence of ftl trimer states [17] is shown 
to lead to an enhanced p-wave scattering cross section in 
the scattering of a 40 K atom and a 40 K- 6 Li dimer. We 
additionally find that the deepest bound trimers have an 
energy spectrum that resembles that of a hydrogen atom 
confined to 2D. 

In the following, we consider a two-component Fermi 
gas confined to quasi-2D by a strong, approximately har- 
monic confinement, Vf,i(z) = ^m^^uj^^z 2 . The gas can 
be considered to be kinematically 2D when the tempera- 
ture and Fermi energy are both smaller than the confining 
frequencies, w-f^ (we set ks = h = 1). The interac- 
tion is characterized by a 3D s-wave scattering length a s 
much larger than the van der Waals range of interatomic 
forces. The two species t, I arc either different hyper- 
fine states of the same atom, or single hyperfine states 
of different atomic species such as 6 Li and 40 K. The low- 
energy scattering of two atoms is described through the 
s-wave scattering amplitude 

fn{q) = ln[l/( g a 2D )]+z^/2' (1) 

where q is the relative momentum. We see here that the 
interparticle interaction is fundamentally different from 
the 3D case: it is energy dependent even at low energies 
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FIG. 1: (a) Illustration of the Skorniakov-Ter-Martirosian 
equation. Wavy (straight) lines are the T-matrix (atom prop- 
agator), respectively, (b) The three-body recombination pro- 
cess. 



and there always exists a bound state with binding energy 
Eb = l/(2jUa2 D ), where the reduced mass = m^ 1 + 
m^ 1 . Equivalently, we can define the scattering T-matrix 
which describes the repeated interparticle scattering at 
total energy E and momentum q, 
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In [{E - q 2 /2(m t + mj,) + iO)/e B ] + in' 



Assuming that the confinement frequencies are identi- 
cal for the two species, then in the limit eb <C uj z , the 
parameter is related to the 3D scattering length by 
«2D = l zy ^rjBe^p{'~^r r j2l z /a s ), with 5 0.905 [2, 20] 
and confinement length l z = l/y / 2fiuj z . For simplicity, in 
this work we restrict ourselves to the 2D limit. However, 
deviations from the 2D limit may be captured by the 
use of a two-channel model where l z plays the role of an 
effective range in the two-body T-matrix [21]. 

We now consider the scattering of a spin-f atom with 
an jj, dimer. The process is illustrated in Fig. 1(a) and 
is described by the Skorniakov-Ter-Martirosian (STM) 
integral equation [22], first developed in the context of 
neutron-deuteron scattering. We let the incoming [out- 
going] atom and dimer have four-momenta (k, e^t) an d 
(-k, E - e k t) [(P, e P f) and (-Pi E - £ pt)L respectively, 
with total energy E = k 2 /2/i3 — eb to ensure that the 
dimer is on-shell. Here we define eict,4. = ^ 2 /^ m t-l and 
the atom-dimer reduced mass p,^ 1 — (m^+m^) -1 +m^ 1 . 
With these definitions, the STM equation in 2D reads 



fe(k,p) = -£h(h,p) ge(k,p) 



qdq ge{p,g)fe(k,g) 
2tt q 2 -k 2 - iO 



(2) 

The atom-dimer scattering preserves angular momentum 
and we have used this to decouple the scattering ampli- 
tude / into partial waves, denoted by £ — for s-wave, 
£ = 1 for p-wave, etc. The projection of the spin-J, atom 
propagator onto the tth partial wave is given by 



gi{p,q) 



cos(£cj))d(j)/2TT 



(3) 



with tf> the angle between p and q. We define the func- 
tion h(k,p) = (k 2 — p 2 )T(p, E — e p ^) to separate out the 
simple pole of the two-particle propagator occuring at 
|k| = |p|. For completeness we also include the factor 
£ describing the quantum statistics. It takes the value 
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FIG. 2: (Color online) s- and p-wave elastic scattering cross 
sections for the fermionic fti system for three different mass 
ratios. Symbols correspond to the low-energy asymptotic be- 
havior, Eq. (5). Note the log-scale for the cross section in (c). 



— 1 when the spin-f particles are fermions, +2 for three 
identical bosons, and +1 in the case of distinguishable 
bosons with negligible f _ t interaction. For a derivation 
of the STM equation in 2D, the reader is referred to the 
Supplemental Material. 

We first discuss the elastic atom-dimer scattering prop- 
erties. The elastic scattering cross section in partial wave 
£ takes the form [23, 24] 
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Note that this has dimensions of length in 2D, and Se : o 
is a Kronecker delta [43]. The partial wave scattering 
amplitudes are related to the solutions of Eq. (2) by 
fe(k) = fi(k,k). Figure 2 shows the s- and p-wave 
scattering cross sections as a function of collision energy, 
-E'coil = fc 2 /2/i3, for equal masses and for the mass ratios 
corresponding to a 6 Li- 40 K heteronuclear Fermi mixture. 
We see that s-wave scattering always dominates when 
to^ < mj,, while p-wave scattering becomes dominant for 
the K-KLi atom-dimer scattering when E co \\/eb > 0.2. 
This is related to a p-wave trimer state crossing the atom- 
dimer continuum at the mass ratio 3.34, leading to a res- 
onant enhancement of p-wave scattering. A similar en- 
hancement of p-wave scattering was predicted in 3D [25] 
for the K-KLi system due to the appearance of a trimer 
at a mass ratio of 8.2 [26]. The presence of trimers in 
higher odd partial waves in 2D at large mass ratio [17] 
will likewise lead to a resonantly enhanced atom-dimer 
interaction. 

In the low-energy limit where fca2D 1, the atom- 
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dimer s- and p-wave amplitudes may be expanded as 
2tt 



/.(*) 



ln[l/(fca ad )]+i7r/2' 



f P (k) « l Sad k\ (5) 



with the atom-dimer scattering length and surface, a a d 
and Sad, describing the s- and p-wave low-energy behav- 
ior, respectively. Matching in the regime fc<22D <C 1, we 
find 
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2.92 
= { 1.54 
0.44 



mf/mi = 1 
rrif/mj. = m K /m u 
mf/mi = m u /m K 



The asymptotic behavior is depicted in Fig. 2. A sur- 
prising feature is that the s-wave scattering cross sec- 
tion is completely determined by a a d for m-j- < toj., 
even though all coefficients in the low-energy expansion 
of the scattering amplitude would be expected to have 
a characteristic scale of <Z2D- A similar result appears 
in 3D for = [27]. The atom-dimer scattering 
length for equal masses has also been estimated from a 
QMC calculation of the pairing gap in the BEC regime, 
but they instead obtain a a d ~ l.7a 2 r> [28]. This dis- 
agreement with our exact few-body calculation is most 
likely because the equation of state has a complicated 
dependence on a a d, making the determination of a a d 
difficult. In the limit of large mass ratios, we find 
fl ad ~ 0.845a2D m(-\/2e 7 / 2 m^//i), with the Euler constant 
7 w 0.577. A logarithmic dependence on mass ratio was 
also found in the 3D case [29]. 

We now turn to trimers, which formally appear as an 
energy pole in Eq. (2) for odd partial waves. The spec- 
trum of p-wave trimers as a function of mass ratio is 
shown in Fig. 3. The appearance of trimers at large mass 
ratios in 2D may be elucidated by applying the Born- 
Oppenheimer approximation [30], which shows that the 
effective potential between heavy fermionic atoms medi- 
ated by the light atom at distances < a2D is similar to the 
electron potential in a hydrogen atom confined to 2D. To 
see this, assume that the state of the light atom at posi- 
tion r adiabatically adjusts itself to the positions ±R/2 
of the heavy atoms. Atom-dimer scattering in odd partial 
wave channels is described by the symmetric light-atom 
wavefunction [29] 

^ R (r) cx K (k(R) ■ |r - R/2|) + K (k(R) ■ |r + R/2|), 

where the modified Bessel function of the second kind 
K (nr) is the decaying solution of the free Schrodinger 
equation with energy e(R) = — /t(i?) 2 /2m^. The sin- 
gularities at the positions of the heavy atoms satisfy 
the Bethe-Peierls boundary condition in 2D: For f = 
r±R/2 this is [r^Yf Mr^o = V ln(fe 7 /2a 2D )- Using 



the asymptotic form, K (x) 



ln(xe 7 /2), we then 



obtain the condition In (— ^p-J = 2/<o 



e(R) R 




FIG. 3: (Color online) Trimer energies (solid lines) as a func- 
tion of the square-root of the mass ratio. Our results match 
the p-wave trimer energies in Ref. [17]. Dashed lines are the 
hydrogen-like spectrum, Eq. (6). 



In the second stage of the Born-Oppenheimer approx- 
imation, the Schrodinger equation of the heavy parti- 
cles is solved using e(R) as the effective interaction po- 
tential. In the limit R <C a2Dj the potential becomes 
e(R) Ri Thus, the spectrum of the deepest 

bound trimer states is hydrogen-like and given by the 
well-known result (appropriately shifted by the dimer 
binding energy) [31] 
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3^m l 2(n+ 1/2) 2 



SB, 



(6) 



with integer quantum number n > £. This also implies 
that deeply bound trimers with different £ are degenerate, 
as was found in Ref. [17]. 

A similar scenario for large mass ratios was recently 
predicted [32] in the context of a 3D Fermi gas close to a 
narrow interspecies Feshbach resonance, characterized by 
a large effective range R* . The effective interaction be- 
tween the two heavy fermions mediated by the light atom 
goes like — 1/R 2 in the range R* <C R <C a s , while at 
shorter ranges, R <C R* , this behavior is replaced by an 
attractive 1/R potential [33] , leading to a crossover from 
an Efimovian spectrum for the weakly bound trimers to 
a hydrogen- like spectrum for the deepest trimers [32]. 

Referring to Fig. 3, for the deepest bound states, we 
find very good agreement with the hydrogen-like spec- 
trum (6). At a given mass ratio, the potential only sup- 
ports a finite number of bound states proportional to the 
number of nodes of the heavy-atom wavefunction that 
fit in the hydrogen-like part of the potential, R < oi2D- 
We estimate this number by noting that in this regime 
the wavefunction of heavy atoms is proportional to the 
Bessel function J 2 e{2\/e~~ / (m^/m^)R/a2n)- Since the 
wavefunction acquires an additional node each time the 
argument increases by 7r, the number of trimers is pro- 
portional to \J m-j- / 77i^ . This feature is clearly observed 
in Fig. 3. 

Turning now to inelastic scattering, three-body recom- 
bination is the process whereby two atoms bind into a 
shallow dimer with the released energy carried away by a 
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FIG. 4: (Color online) The recombination rate as a function of 
energy for (a) m-f = and (b) the Li-K heteronuclear mix- 
ture. In the low energy limit E/eb <C 1, the result for identi- 
cal bosons matches the purely s-wave calculation in Ref. [18] . 



third atom. To extract the recombination rate from the 
atom-dimer scattering amplitude, we consider the recom- 
bination process in reverse. This is illustrated in Fig. 1(b) 
which shows how the recombination process is described 
by the same diagram as inelastic atom-dimer scattering - 
a process in which no dimer remains after the scattering. 
The total recombination rate is therefore proportional 
to the atom-dimer inelastic scattering cross section [34] 
K{E) = v^^-J^i^T^i^)^ where the atom-dimer rel- 
ative speed is w a d = k/p,^ and the sum is over all partial 
waves. The phase space of an atom-dimer pair 4> a( j and 
three atoms $ aaa are given in the Supplemental Material. 
The inelastic scattering cross section er™ 01 is obtained by 
subtracting the elastic from the total cross section, which 
is related to the imaginary part of the scattering ampli- 
tude through the optical theorem [23, 24] 
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af t (E) = -%[ft(k)](2-5 i!0 ). 



Thus we arrive at the recombination rate 
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(8) 



with the degeneracy factor p = 2 if f^-l, i.e. for fermions 
or heteronuclear bosons, and p = 3 for identical bosons. 

Figure 4 shows the energy-dependent recombination 
rate for equal masses and for the Li-K fermionic mixture. 
We see that in the mass-balanced case, the fermionic re- 
combination rate peaks at E ps 3eb instead of E ps 8eb 
as in the case of identical bosons. Moreover, we find that 
the fermionic recombination rate is bounded by the rates 
for identical and heteronuclear bosons, with the fermionic 
rate being comparable to that of identical bosons in the 




FIG. 5: (Color online) The event rate constant a for fermions 
with m-f — mi as a function of (inverse) 3D scattering length. 



low energy limit. This is manifestly different from the be- 
havior in 3D, where the low-energy recombination rate 
approaches a constant in the three-boson problem [35] 
while it is suppressed in Fermi systems due to Pauli 
blocking with K(E — !• 0) ~ E [29]. Surprisingly, we 
find that the ei-wave channel dominates recombination 
of K-K-Li fermionic atoms and of three identical bosons 
for most of the plotted energy range, whereas the p-wave 
channel is dominant for equal mass and Li-Li-K fermions, 
as well as for heteronuclear bosons. In all cases we find 
that K(E) approaches zero faster than E, which is a 
purely two-dimensional feature. 

To obtain an experimentally accessible quantity, we 
average K{E) over energy with an appropriate thermal 
distribution. This yields the event rate constant per unit 
volume a, defined such that the number of recombina- 
tion events per unit volume per unit time is anin± for 
the process t + t + + Tl i n a two-component Fermi 
mixture. When the kinetic energy of the resulting parti- 
cles is greater than the height of the trap, the loss rate 
is given by = — 2cv^n|n^ — a^njn^. At sufficiently 
high temperatures, we employ the Boltzmann distribu- 
tion which gives [18, 36] 



a(T) 



r °° dEEe- E ' k * T K{E) 
p\j™dEEe- E / k BT ■ 



(9) 



This quantity is illustrated in Fig. 5 for a range of typical 
temperatures [44]. We see that for increasing tempera- 
tures, the peak value decreases as a pca k ~ T _1 and the 
peak position, [/ z /a s ] pca k ~ — hi(T), shifts from the BCS 
side of the 3D resonance towards unitarity. 

At present, heteronuclear fermionic Li-K mixtures [37, 
38] are promising candidates for the observation of res- 
onantly enhanced atom-dimer scattering and trimer for- 
mation. Avenues towards observing the hydrogen-like 
spectrum predicted at large mass ratio include using the 
recently predicted Feshbach resonance in fermionic Li- 
Yb mixtures [39] as well as the application of a species 
selective optical lattice to enhance the effective mass ra- 
tio [40]. 

We also expect our results to be important for ongo- 
ing experiments on quasi-2D Fermi gases. In particular, 
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three-body recombination will limit the stability of re- 
pulsive Fermi gases and potentially generate three-body 
correlations, thus impacting any experiment seeking to 
realize itinerant ferromagnetism. Moreover, it would be 
interesting to investigate how three-body losses might 
evolve into two-body losses in the presence of a Fermi 
sea [41, 42]. 
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Skorniakov— Ter-Martirosian equation two dimensions 

In 2D, the Skorniakov-Ter-Martirosian integral equation is given by (see Fig. lb in the main text), 



T(k,/c ;p,p ) =C 



ZG i {--k-p,E-k -po) 
d 2 q dq 



+ i 



(2tt) 2 2tt 



G t (q, q )Gi{-p -q,E-p - q )T(-q, E - q )T(k, k ; q, q ) 



(1) 



Here, T(k, fc ;p,po) is the atom-dimer scattering T- matrix. We let the incoming [outgoing] atom and dimer have 
four-momenta (k, e^) and (— k, E — e^) [(p, e p -f) and (— p, E — e p -j-)], respectively. The free atom propagators and 
the two-particle T-matrix are given by 



G t ,i(k, fc ) = 
T(q,<?o) = 



1 

fc - fc 2 /2m t ^ + i0' 

2tt/^ 



- In [(q - q 2 /2(m 1 - + mi) + iO)/e B ] + in ' 



(2) 
(3) 



with /! 1 



m | . The residue at the pole of the two-particle T-matrix is used to extract the dimer wavefunction 



renormalisation Z = 2t:eb/^- The factor £ describes the quantum statistics and is equal to —1 for in the fermionic 
case, +1 for heteronuclear bosons, and +2 for identical bosons. 

Integration over q can be done by closing the contour in the lower half plane. This encloses one pole coming from 
G-f(q, q ) at q a = e q f. Next, we take the partial wave projection and use the on-shell conditions: fc = e^, po = e p -j- 
and E = k 2 /2fi 3 — eb with ^3 1 = m^ 1 + (m-j- + m.).) -1 . The equation (1) becomes 



Te(k; P ) = ( 



Zge(k,p) + J ^-ge(p,q)T(q,E - e^)T e (k; q) 



The partial wave projection in 2D is defined as follows 



Ti(p,q) ■■ 

where <f> is the angle between p and q. Similarly, 

^ d<f> 



2tt 



cos(^)T(p,q), 



9e(p 



r 

Jo 



2tt 



cos{£(f>)Gi{p - q, E - e pt - e qT ) 



= I Sign(E-p 2 /2^-q 2 /2^ y/(E - p 2 /2^ - q 2 /2^ + 1O) 2 - (pq/mj- 



{ [Y{p, q)X e (p, q) - X t -i(p, q)] g (p, q) - (m l /pq)X e (p, q) 
where we define, for any non-negative integers I 



if £ = 
if £> 1 



(4) 

(5) 
(6) 



Xi(p,q) 



(Y(p,q) + y/Ytfaq)-!)' - (Y(p,q) - ^Y 2 (p, q) - l) ' 



^Y 2 {p,q) 



Y( P ,q) = ^(E-l±^ 
pq \ 2fi 



(7) 
(8) 
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The elastic atom-dimer scattering amplitude from incoming momenta ±k to outgoing ±k' is given by /(k, k') = 
2/i 3 T(k, Ck^; k', ekf) with |k| = |k'|. Then we define functions / and h: 



ft(k,q) _T(q,E-e vt ) 



<r iO z T ^ 

h(k,q) = (k 2 -q 2 )T(q,E-e clt ), 



(9) 
(10) 



such that the partial atom-dimer amplitudes are given by fe(k) = fe{k, k). The functions ft obey a set of decoupled 
equations for each partial wave: 



fe(k,p) = -(h(k,p) ge(k,p) + J 
The scattering amplitude is related to the phase shift by 



qdq ge(p,q)fe(k,q) 



2ir q 2 -k 2 - iO 



ft(k) = 



cot(Si) + i ' 



with the low-energy expansion 



mtf^-J (2A)Ml/(fcaad)] if £ = 



(11) 



(12) 



(13) 



valid for ka2B <C 1- 



Phase space 

The TV-particle phase space integral is defined such that 



where indistinguishability of particles gives rise to the factor r\ < 1. Using this definition, the atom-dimer and 
three-atom phase space is given by 

$ad = A*3/27T, ^aaa^-T^ " " • (14) 

For three particles, rj = 1/ p\ where p is the number of identical particles. The phase space factors were derived in the 
bosonic case for equal masses in Ref. [18]. 



